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The dynamics of individual colloidal particles in random potential energy landscapes were investigated experimentally and by 
Monte Carlo simulations. The value of the potential at each point in the two-dimensional energy landscape follows a Gaussian 
distribution. The width of the distribution, and hence the degree of roughness of the energy landscape, was varied and its effect on 
the particle dynamics studied. In the experiments, the energy landscapes were generated optically using a holographic set-up with 
a spatial light modulator, and the particle trajectories followed by video microscopy. The dynamics are characterized using, e.g., 
the time-dependent diffusion coefficient, the mean squared displacement, the van Hove function and the non-Gaussian parameter. 
The dynamics are initially diffusive, show an extended sub-diffusive regime at intermediate times before diffusive motion is 
recovered at very long times. Compared to the dynamics in a one-dimensional potential energy landscape, the localization at 
intermediate times is weaker and the diffusive regime at long times reached earlier, which is due to the possibility to avoid local 
maxima in two-dimensional energy landscapes. 
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1 Introduction 

The Brownian motion of colloidal part icles is one of the clas- 
sical phenomena in statistical physics. I^lEI in real situations, 
the particles often do not move freely , but their dynamics are 
modified by an external potential. Especially a random po- 
tential, and thus Brownian motion in the presence of disor- 
der, leads to interesting transport phenomena. ^ ^ Up to now, 
the dynamics in random potentials have been studied mainly 
by theory and computer simulations.^ lO ii 12 13 14 15 Theoret- 
ical models include the random barrier model,EI the random 
trap model, the random walk with barriers and the con- 
tinuous time random walk^^ as well as studies of diffusion 
in a rough potential and in materials with defects like zeo- 
lites. In particular the long-time limit has been investigated 
for different realizations of random potentials.^ ^ In contrast, 
less is known on the intermediate regime and the time needed 
to reach the long-time limit. To our knowledge, only very 
few systematic experimental tests o f theore tical and simula- 
tion predictions have been performed.E^E^lSlNevertheless, the 
theoretical predictions have been applied successfully to ex- 
perimental data and the concept of particles diffusing through 
an energy landscape has proven very useful in understanding 
very different phenomena. This includes particle diffusion 
in inhomogeneous medi a (e.g. single molecule dynamics in 
porousgelsEH or in cellsEEES), the dynamics on rough sur- 
f aces, 1232^ the dynamics of particles moving along the walls 
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between magnetic domains, the dynamics of independent 
charge carriers in a conductor with impurities (in the param- 
eter ran ge wh ere conduction can be modeled as a classical 
process).E^^Furthermore, some processes can be represented 
by diffusion in the systems' configura tion space, fo r example 
vitrification leading to glassy systems 3 Q|3i|32 | 33 | 34|35 | pj-Qtein 
folding. Often diffusion in a potential energy land- 

scape represents a crude approximation only, but it can nev- 
ertheless provide a useful first description of the effect of dis- 
order on the dynamics. Disorder may modify the value of the 
diffusion coefficient or it may alter Brownian motion leading 
to anomalous diffusion. Which effect dominates depends not 
only on the specific process, but also on the time scale of in- 
terest. 

An external potential can be imposed on a polarizabl e col- 
loidal particle by exposing it to a light field.E^E^^^SI Light 
exerts different forces on particles, if their refractive index 
differs from (typically exceeds) that of the solvent: a scatter- 
ing force or 'radiation pressure', which pushes particles along 
the laser beam, and a gradient force, which a ttracts particles 
toward regions of high light intensity.S^E^Sl a classical ap- 
plication of this effect are optical tweezers which are used 
to trap individ ual particles by a tightly focused laser beam. 
| 4i|42 | 43 | 45 | 46|47 | Furthermore, above a certain light intensity, a 
periodic light field can induce a disorder-order transition in 
a two-dimensional charged colloidal system, known as light- 
induced freezing, and, if the intensity is increased further, the 
induced crystal melts in to a modulated liquid, called light- 
induced melting .I^HEmni in addition to the particle arrange- 
ment, the particle dynamics can be affected by periodic and 
random light fields, resulting in anomalous diffusion. Light 
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fields hence provide a means to manipulate the spatial arrange- 
ment and dynamics of colloidal particles. 

Recently, we experimentally realized one-dimensional ran- 
dom energy landscapes and periodic potentials'^"^ using 
laser light fields and studied the dynamics of individual parti- 
cles in these potentials. Here, this is extended to the dynam- 
ics of individual colloidal particles in two-dimensional ran- 
dom potentials. In our experiments and simulations, the val- 
ues of the two-dimensional random potential were drawn from 
a Gaussian distribution, whose width £ represents the degree 
of roughness of the potential and in the experiments was con- 
trolled by the laser power P. The static properties of the po- 
tential were determined quantitatively. Furthermore, the tra- 
jectories of individual p articles i n this potential were followed 
using video micro scopy,EElEl and compared to our simula- 
tion results. The dynamics were characterized by, e.g., the 
time-dependent diffusion coefficient, the mean squared dis- 
placement (MSD), the non-Gaussian parameter, and the van 
Hove function. The dynamics are initially diffusive but then, 
at intermediate times, show an extended subdiffusive regime 
before diffusive behaviour is reestablished at very long times. 
Our findings are compared to the particle dynamics in one- 
dimensional random potentials and periodic potentials. 
EH Two dimensions allows particles to bypass large barriers. 
Therefore, the particle dynamics are controlled by minima and 
saddle points instead of minima and maxima. Moreover, com- 
pared to periodic potentials the barriers have different heights. 
This significantly affects the particle dynamics, for example 
the dependence of the dynamics on the degree of roughness of 
the potential is modified. 

2 Materials and methods 
2.1 Sample preparation 

Each sample consisted of surfactant-free sulfonated 
polystyrene (PS) particles with a radius R = lA \xm and 
polydispersity 3.2 % (Interfacial Dynamics Microspheres 
& Nanospheres) suspended in heavy water (D2O), so that 
the particles cream rather than sediment. Stock solutions 
of the particles were diluted to result in an area fraction of 
the creamed sample, a < 0.10, which represents a compro- 
mise between negligible particle-particle interactions and 
reasonable statistics. Area fractions were estimated from 
micrographs according to cr = 7lR^Np/A with A^p and A being 
the number of particles and the area covered by the light field, 
respectively. 

The heavy water (D2O) was de-ionised by stirring with ion 
exchange resin to increase the particle-glass repulsion and thus 
reduce the fraction of particles sticking to the glass surface. 
To further reduce sticking, all glassware was sonicated in 2% 
Helmanex II solution at about 60 °C and then rinsed with Mil- 



lipore water and dried in air prior to use. Each sample cell was 
constructed from a microscope slide and three cover slips, two 
used as spacers (number with thickness 0.085 — 0.13 mm, 
supplied by VWR) with a gap between them and the third 
on top to create a narrow capillary (number 1 with thickness 
0.13 — 0.16 mm, supplied by VWR).^^ Thin cover glasses 
were used as spacers to allow imaging of the creamed parti- 
cles using a high resolution objective with a working distance 
of 0.13 mm. The sample chamber was filled using capillary 
action and subsequently sealed with UV glue. 



2.2 Light field generation 

The set-up contains a laser with a wavelength of 532 nm (Ven- 
tus 532-1500, Laser Quantum). Its beam is expanded and then 
reflected from a spatial light modulator (Holoeye 2500-LCR). 
Subsequently it is directed through two telescopes to reduce 
its diameter and reflected off three mirrors to steer it through 
an inverted microscope (Nikon Eclipse 2000-U) into the sam- 
ple ii9|45, One of the mirrors is a dichroic mirror to introduce 
the beam into the microscope beam path and to use the mi- 
croscope objective (60 x oil immersion, numerical aperture 
NA 1.4, Nikon) to image the light field into the sample plane. 
The beam passes upwards through the sample and hence, due 
to radiation pressure, pushes the particles against the top of the 
cell, which reinforces the creaming of the particles. A notch 
filter in the imaging path prevents laser light from reaching 
the ocular or camera. To aid alignment, the notch filter can be 
removed and the sample replaced by a mirror, so that the light 
intensity distribution in the sample plane can be imaged using 
the microscope. 

A kinoform (phase hologram) was calculated using the 
Gerchberg-Saxton iterative algorithm (Fig. [T]\) and dis- 
played in the centre of the spatial light modulator. The ki- 
noform corresponds to a homogeneous disc surrounded by a 
ring to prevent particle movements into and out of the disc. 
The Fourier transform of the kinoform is, as expected, a ho- 
mogeneous disc surrounded by a ring (Fig. [Tp). In order to 
account for the angle at which the laser impinges on the spa- 
tial light modulator (22.5°), the disc and rin g are a factor 
1/ cos (22.5°) = 1.08 taller than they are wide.ESEl xhe ob- 
served light field intensity I{x,y) (Fig.[Tp) corresponds to the 
disc of the Fourier transformed kinoform. Indeed, the illumi- 
nation is overall flat but, crucially, has some fluctuations due 
to the finite size and pixelation of the light modulator. These 
fluctuations are exploited in the following. Furthermore, there 
is a bright 0* -order peak in the centre. Using this peak, a par- 
ticle was trapped and used to monitor any drift of the set-up. 

Global drifts were found to be negligible during individual 
measurements (up to 4 h). 
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Fig. 1 (A) Kinoform calculated by applying the Gerchberg-Saxton 
algorithm to a homogeneous disc surrounded by a ring and (B) its 
Fourier transform. (C) Micrograph of the observed intensity I{x,y) 
of the disc taken at very low laser power P < 0.2 mW. (D) Potential 
U (x,_y) as experienced by a point-like test particle obtained by 
convoluting I{x,y) with the volume of a spherical particle with 
radius R = \ A \im =12.7 pixels. 

2.3 Video microscopy and particle tracking 

The samples were observed using the inverted microscope. 
Micrographs were recorded using a CMOS camera (PL- 
B742F, Pixelink). Particle coordinates were extracted from the 
time series of micrographs and the trajectories determined us- 
ing IDL routines. Typical measurement times were 2 to 4 h. 
Particles which were stuck to the glass were identified by com- 
paring the particles' short-time friction coefficient ^j, i.e. the 
inverse mobility, determined from the mean squared displace- 
ment, to the expected bulk value = GTirfR, with the solvent 
viscosity r] = 1.19 x 10~^ Pas at room temperature. Particles 
with > 20(^0 were declared stuck and removed from the 
analysis. Typically, one particle was stuck to the glass in the 
field of view, which contained about 20 particles. 

For identical conditions, measurements at different posi- 
tions in the sample yielded very similar results, despite slightly 
different particle area fractions a. This reproducibility al- 
lowed us to average several independent measurements to im- 
prove statistics. 

2.4 Monte Carlo simulations 

The Monte Carlo simulations were performed on a 4096 x 
4096 square lattice with the lattice points separated by in 



Fig. 2 Some region of the spatially correlated Gaussian potential 
energy landscape U (x, j), obtained by convolution of a spatially 
uncorrected Gaussian energy landscape with the particle volume. It 
is used in the Monte Carlo simulations and reflects the potential felt 
by a particle (Fig.|Tp). 

both directions, where we have set = 1. The potential val- 
ues at the lattice points, U (x, j), were produced using a Box- 
Muller algorithm generating numbers which are Gaussian dis- 
tributed with zero mean and standard deviation £. The poten- 
tial U (x, was convoluted with the particle volume to obtain 
the potential U {x^y) felt by a point-like test particle 

Y,Y,U{x—kAs,y—lAs) a{k,l) 
U{x,y) = ^ , (1) 

where the double sum runs over the projected particle, i.e. + 
< with kAs and I As the distances from the particle centre 
in the two directions and R = mAs the radius of the particle. 
The volume of the particle is represented by 

a{k, I) = 2^J{m^-k^-fi) . (2) 

As a compromise between negligible discretization effects and 
viable computation time, we have chosen m = 32 and thus 
-32 <k,l< 32 with k^^f <n?. 

The convolution leads to a potential U (x^y) (Fig. [2]), which 
is smoother than U (x, j). Its values follow the same Gaussian 
distribution, albeit with a spatial correlation decaying on the 
length scale of the particle size. It is supposed to resemble the 
potential energy landscape experienced by a colloidal particle 
in the light field (Sec. [31]). 

Once the potential energy landscape U{x^y) was fixed, a 
particle was set on a randomly chosen lattice point. During the 
simulation, a direction is chosen randomly and, depending on 
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the energy difference AC/ to the neighbouring lattice point, the 
particle is moved in any case if AC/ < 0, or moved with a finite 
probability exp (— AC//^b^) if AC/ > 0. By averaging over 
1024 different positions of the particle, representative aver- 
ages can be determined. For each Monte Carlo run, the short- 
time diffusion coefficient Dq and the related Brownian time 
= R^/4Do were calculated. In analogy to the experiment, 
data were acquired up to 10^ ^b- This yielded particle trajecto- 
ries as in the experiments. Thus, the different parameters, such 
as the mean squared displacement, were determined as in the 
experiments and are described below. It turned out that within 
statistical errors the results for different realizations of the po- 
tential energy landscape U (x^y) are identical. Moreover, as in 
the experiments, separate simulations were performed for dif- 
ferent values of the degree of roughness Ok^T <£ <5kBT in 
steps of 0.25 k^T to investigate its effect on the dynamics. 

3 Results and discussion 

We studied the behaviour of individual colloidal particles 
in two-dimensional random potential energy landscapes. At 
first, the properties of the experimentally created energy land- 
scapes are presented. Then, the particle dynamics in these 
energy landscapes are discussed and compared to the results 
of our Monte Carlo simulations. Finally, our experimental 
and simulation results are contrasted to the dynamics in one- 
dimensional random and periodic potentials. 

3.1 Properties of the optically generated random poten- 
tial 

A realization of the light field at very low laser power is dis- 
played in Fig lip . The light field interacts with polarizable 
particles. '^I'^^Eirxhe polarizable particle volume is taken into 
account by convolving the local light intensity I{x^y) with the 
particle volume. The effect of the light field on the particle is 
then represented by an external potential U{x^y) as felt by a 
point-like test particle (Fig.[Tp). Due to the convolution with 
the particle volume, the potential is not only proportional to 
the laser power P, but also to the particle volume, U oc PR^ . 

To determine the characteristic length scales of the light 
field intensity I{x^y) and of the potential felt by a point-like 
test particle U {x^y), the spatial correlation functions were de- 
termined. Based on the two-dimensional intensity distribution 
I{x,y), the azimuthal average of the two-dimensional spatial 
correlation function {^{x^y))^ were calculated. The spatial 
correlation of the light field intensity I{x^y) decays on a short 
length scale compared to the particle size. However, the con- 
volution with the particle volume introduces a length scale, 
namely the particle diameter 2R. The spatial correlation of 
the potential U (x, y) , which was similarly determined, indeed 
decays on a characteristic length of 2R (Fig. [3]). 
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Fig. 3 Azimuthally averaged spatial correlation functions of the 
laser intensity I{x,y) (dashed line) and of the potential energy 
landscape felt by a point- like test particle U {x,y) (solid line) as 
determined from Fig.[Tp and D, respectively. 

Based on the observed light intensity I{x^y) and potential 
energy landscape U{x^y) (Fig.[Tp,D), the distributions of the 
light intensity values p{I) and potential values p{U) were de- 
termined (Fig.|4]). The distribution p{I) follows a Gamma dis- 
tributional 
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where / > and r(^) is the Gamma function. A fit to the ex- 
perimental p{I) yielded /: = 3.1 ±0.1 (Fig.|4|\), correspond- 
ing to a 3D speckle pattern.*^^*^ The distribution p{U) can be 
described by a Gaussian distribution 
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with the average (C/) and width or standard deviation £ 
(Fig. |4^). Due to the convolution with the particle volume, 
U{x^y) represents a weighted average of several independent 
(random) values of I{x,y) and thus p{U) has a significantly re- 
duced width compared to p{I). The width £ characterizes the 
degree of roughness of the random potential C/(x,j), which 
is controlled by the laser power P, but cannot be easily de- 
termined experimentally. Thus, to establish a quantitative re- 
lation between the roughness £, used in the simulations, and 
the laser power P, applied in experiments, the experimental 
potential energy landscape needs to be calibrated. This was 
achieved by a direct comparison of the experimental and sim- 
ulation results, namely of the time-dependent diffusion coeffi- 
cient D{t) at very short and long times, and will be described 
below (Sec. |3.3| ). The calibration resulted in an approximately 
linear relation between £ and P, which might saturate for large 
P(Fig.|5j. 
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Fig. 4 Distribution of (A) values of the intensity of the Ught field, 
p{I), and (B) values of the potential as felt by a point-like test 
particle, p{U) based on the observed intensity I{x,y) and potential 
U (x,_y) shown in Fig.[Tp and D, respectively. Red lines are fits 
based on a Gamma and Gaussian distribution, respectively. 



3.2 Dynamics in the random potential - experiments 

The effect of two-dimensional random energy landscapes on 
the particle dynamics is qualitatively illustrated in Fig. [6] Out- 
side the light field (white background), particles undergo free 
diffusion, exploring a large area. Separated from this region 
by a large barrier (white/green rings) is the two-dimensional 
random light field (green disc). Within the random potential, 
the excursions of the particles are limited and the particle dy- 
namics hence are slowed down. The particles remain longer 
at some positions, which correspond to local minima of the 
potential. For a potential with a larger degree of roughness e, 
i.e. a larger width of p{U), this effect is more pronounced with 
particles being more efficiently trapped and hence exploring a 
smaller region. 

Based on the particle trajectories, different statistical prop- 
erties were computed to characterize the particle dynamics. 



Fig. 5 Standard deviation e of the distribution of potential energy 
values, p{U), as a function of laser power P. 
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Fig. 6 Trajectories of particles undergoing diffusion in a 
two-dimensional plane, part of which contains a random potential 
(green background) which is separated by a barrier (white/green 
rings) from the surroundings (white background). Particle radius 
R= \ A |am, particle surface fraction a = 0.04, laser power 
P = 1.32 W corresponding to a standard deviation e = 2.5 k^T, and 
a recording time T = 3.8 h. Coordinates are given in fim. 



We found identical behaviour along the x- and j-directions 
as expected for an isotropic system. The dynamical prop- 
erties were hence determined as a function of the distance, 
Ar = [(Ax)^ + (Aj)^]^/^, where distances are scaled by the 
particle radius R = lA jam and times by the Brownian time 
tB =R^/ (4Z)o) = (6.3 ± 0.1) s (values refer to the experimen- 
tal system; Dq has been determined in the absence of a ran- 
dom potential, i.e. e = 0, but in the vicinity of the water-glass 
interface). This renders the data independent of the specific 
experimental conditions. Moreover, the statistical properties 
were obtained by averaging over different particles, which are 
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Fig. 7 Particle residence time distribution ^/,£(0 representing the 
probability that it takes a particle a time t /t^ to travel at least a 
distance / in a random potential with standard deviation e. All 
curves are smoothed by a moving five-points average. (A) ^[^^{t) 
for different length l/R (as indicated) and e = 2.5 k^T, scaled plot 
as inset. (B) ^/^^(O for l/R = l and different e (as indicated). 



non-interacting and thus independent, and waiting times to 
improve statistics. Since initially the particles were homoge- 
neously distributed, in the experiments and simulations, but 
tended toward a Boltzmann distribution in the course of the 
experiment or simulation, the average over waiting times de- 
pends on the total measurement or simulation time T, which 
was r 1000 ^B. 

Depending on the particle positions, the particles experi- 
ence various potential values U {x^y) and are trapped for differ- 
ent times, reflecting the different heights of the saddle points 
to the neighbouring minima. The time t required to explore 
at least a distance / in a potential with roughness £ has been 
determined and the particle residence time distribution ^/,e(0 
determined. To explore a distance / by free diffusion with 
diffusion coefficient Dq, on average the time = ^^/(4/)o) 
is required. To explore larger distances / and/or in the pres- 
ence of a random potential, on average larger times are re- 
quired (Fig.[7]A). For short distances / < 2R, i.e. within a min- 
imum, ^/,e(0 does not significantly depend on e. However, to 
travel a characteristic length 2R, which corresponds to the typ- 
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Fig. 8 Distribution of particle displacements Ar within time t, 
P{Ar,t) for (A) free diffusion, i.e. without potential (e = 0), scaled 
P{Ar,t) as an inset, (B) dynamics in the presence of a random 
potential with roughness £ = 2.5 k^T for different times t (as 
indicated) and (C) with different roughnesses £ (as indicated) for 
time ^ = 50^B- 

ical minimum-minimum separation (Fig.[3|, on average a time 
t = 4tB = R^/Dq is required and, in the presence of a random 
potential, the distribution ^/,e(0 depends on the roughness £ 
(Fig. [7^), since a barrier or saddle point, whose average height 
depends on £, has typically to be crossed to travel a distance 
1>2R. 

The probability distribution of particle displacements Ar, 
i.e. the self part of the van Hove function, P(Ar, t), at different 
delay times t is calculated based on the trajectories by averag- 
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ing over all waiting times and particles / 

P(Ar,0 = (5(Ar- [r^^o+O - ^(^o)] , (5) 

where r/(^) is the position of particle / at time t. In the case 
of free two-dimensional diffusion, i.e. without any external 
potential, P{Ar,t) follows a Rayleigh distribution, P{Ar,t) ^ 
Ar/{2Dot) exp (— Ar^/4Z)o^), whose width increases linearly 
with time t (Fig. [8}\). In the presence of a random potential, 
P{Arj) changes qualitatively (Fig. [8^). The potential tends to 
trap the particle so that it explores less space and the distribu- 
tions P{Ar^t) get much narrower. This is more pronounced for 
longer times, when the dynamics include barrier crossing. Ac- 
cordingly, at long delay times, the roughness of the potential 
significantly effects P(Ar, t), which becomes narrower with in- 
creasing £ (Fig.[8p). 

The width of the distribution of particle displacements, 
P{Ar,t), can be characterized by the mean squared displace- 
ment (MSD) 

(Ar2(0),^. = ([n(^o+0-n(^o)]')^^. (6) 

where {Ax^(t)) = {x^(t)) - {x{t))^ and {Ay^{t)) is calculated 
correspondingly. The particle dynamics, as reflected in the 
MSDs, exhibit three distinct regimes (Fig.|9]A). Both, at short 
times (t/tB ^0.1) and long times (t/tB ^ 30), the particle dy- 
namics are diffusive. At small t, the diffusive behaviour re- 
flects small excursions within local minima and is thus essen- 
tially independent of the roughness £. (Nevertheless, diffusion 
at short times is reduced compared to free diffusion (e = 0), 
because averages over waiting times are considered and 
laser pressure pushes the particles clos er to the water-glass in- 
terface and thus reduces their mobility,!^ ^^ ! ^^ ! ^ ^ I with only a very 
weak apparent dependence on laser power P > and hence 
e > 0; see Sec. |3.3[ ) For large enough t, hopping between min- 
ima becomes important and constitutes a random walk. Thus 
diffusive behaviour is reestablished at long times, although 
with a strongly reduced diffusion coefficient. At intermedi- 
ate t, the MSDs exhibit an inflection point, which becomes 
increasingly pronounced as £ increases. This subdiffusive be- 
haviour is caused by the particle being trapped in local minima 
for prolonged times before it escapes to a neighbouring min- 
ima. Since there is a wide range of residence times (Fig. |7]), 
reflecting barriers of different heights, the subdiffusive regime 
extends over a broad range of times. 

From the two-dimensional MSD {Ar^{t)), the time- 
dependent diffusion coefficient D{t) can be calculated accord- 
ing to 
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Fig. 9 Particle dynamics in a two-dimensional random potential 
energy landscape with increasing degrees of roughness £ (as 
indicated): (A) normalized mean squared displacement 
{Ar{tf) /R^, (B) normalized diffusion coefficient D{t)/Do, (C) 
exponent ^(^) in the relation (Ar^(^)) oc t^, and (D) non-Gaussian 
parameter 0^2 (0 ^ function of delay time t normalized by the 
Brownian time • For clarity, only every fifth data point is plotted 
as a symbol. Black crosses indicate minima and maxima of ^(^) and 
a2{t), respectively. 



where here the dimension d = 2. The three regimes discussed 
above are also reflected in the normalized time-dependent dif- 
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fusion coefficient D{t)/Do (Fig.|9^). Toward very short times, 
D{t)/Do tends toward one (actually slightly below one due 
to the averaging and the radiation pressure discussed above 



and in Sec. |3.3| ). It strongly decreases at intermediate times to 
reach a much smaller value Doo at long times, where hopping 
between minima dominates and diffusion is reestablished, re- 
flected in the plateau of D{t). The long-time regime is exper- 
imentally accessible and hence the asymptotic diffusion coef- 
ficient Doo can be determined experimentally and will be dis- 



cussed together with the simulation results (Sec. |3.3| ). 

In order to characterize deviations from diffusive behaviour, 
in particular the subdiffusion at intermediate times, the expo- 
nent /i in the relation (Ar^(r)) ~ is considered. For free 
diffusion /i = 1, while /i < 1 in the case of subdiffusion. The 
time-dependence of ii{t) is calculated from the slope of the 
MSD in double-logarithmic representation: 



aiog(<Ar^(0),^J 



dlog{t) 



(9) 



The subdiffusive dynamics at intermediate times cause a min- 
imum in ii{t), which becomes more pronounced with increas- 
ing £, but remains at about the same time (Fig.[9p, crosses). 
In contrast, the diffusive behaviour at short and long times is 
reflected in the trend of tending toward one for the two 
limits. 

While the exponent characterizes deviations from dif- 
fusive behaviour, the non-Gaussian parameter a2{t) quanti- 
fies the deviation of the distribution of particle displacements, 
P{Arj), from a Gaussian distribution (Fig.jsj). It corresponds 
to the first non-Gaussian correction to P{Ar,t).^^ Following a 
previous definition:E2 



a2{t) 



{Ar\t)) 



{1^2/d){Ar^{t))l. 
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(10) 



where (Ar'^it)) = {Ax\t)) ^ {Ay\t)) ^ 2 {Ax^ {t)) (Ay^ {t)) 
and (^Ax^{t)) and (A_y^(^)) are defined in analogy to (Ax^(r)). 
The time-dependence of CC2{t) also shows three different dy- 
namic regimes (Fig. [9p). At very short and very long times, 
when the particle dynamics are diffusive, a2{t) ~ 0, while at 
intermediate times 0^2 (^) develops a peak which becomes more 
pronounced and moves to larger ta with increasing £. This 
reflects the broader distribution of barrier heights and hence 
residence times ^/,e(0 at larger £ (Fig.]?]). 

The intermediate regime is hence characterized by subdif- 
fusive dynamics with a minimum in the exponent il{t) and, 
at later times, a maximum in the non-Gaussian parameter 
0C2{t) (Fig.[9p,D). While the position of the minimum in ii{t), 
i.e. the time hardly depends on e, the maximum in a2{t) 
shifts to larger times ta with increasing £ (Fig.fTO). 




Fig. 10 Characteristic times, namely of the minimum in the 
exponent jU(^), i.e. t^, and the maximum in the non-Gaussian 
parameter a2{t), i.e. ta, as a function of the degree of roughness £ 
of the potential energy landscapes from simulations (open symbols) 
and experiments (filled symbols, corresponding to the crosses in 
Fig.[9p,D). Solid lines are guides to the eye. 




Fig. 11 Normalized diffusion coefficient, D(t)/DQ, as a function of 
delay time t for particles in a two-dimensional random potential 
energy landscape as observed in simulations (black solid lines, 
roughness £ as indicated) and experiments (coloured lines with 
symbols, roughness £ as indicated). To account for the effect of 
radiation pressure on the particle mobility, the experiments are 
scaled with an effective Brownian time t^ and an effective diffusion 
coefficient Dq as explained in the main text. 



3.3 Dynamics in the random potential - simulations 

The experimental findings, namely the time-dependent diffu- 
sion coefficient D{t), will be compared to simulation results. 
The degree of roughness £ of the optically-generated poten- 
tial U{x^y) is tuned via the laser power P, but £ cannot be 



8| 



1^ 



determined experimentally. Thus, the relation between the 
roughness £ and the laser power P needs to be determined by a 
comparison of the simulation and experimental findings. The 
diffusion coefficient D{t) exhibits well-developed plateaus at 
early and late times, = D{ts) and Doo, which depend on 
£ (in the simulations) and P (in the experiments). Therefore 



we use the ratio Doo/D^ determined in simulations (Fig. 12 A) 



to assign an £ to the Doo/D^ which have been determined ex- 
perimentally using different P. This yields a relation between 
roughness and laser power, i.e. £{P) (Fig. [5]), which appears 
linear up to large P, where £ starts to saturate. Also other 
procedures have been followed to determine £{P)', they all re- 
sulted in an apparently linear relation between £ and P and 
slopes of £{P) within 20%. Furthermore, the slope is consis- 
tent with a previous calibration of a one-dimensional random 
potential, when taking the different illuminated areas into ac- 
count.IH 

In addition, the friction coefficient of the particles, ^ , im- 
plicitly depends on P. Due to hydrodynamic effects, the 
friction coefficient varies with the particles' distance from 
the water-glass interface. ^^^^^^ The distance is controlled 
by a bala nce between the repulsive particl e-wall interaction 
| 63|64|65 | 66| radiation pressure,'^^!!!!^] pushes the 

particles toward the glass slide and depends on the laser power 
P. Thus, a finite P > will lead to an increased friction 
coefficient and hence reduced diffusion coefficient. At short 
times the diffusion coefficient tends to a value, D^, which can 
hence be used to guide the correction. Although the dynam- 
ics at short times are hardly affected by the random poten- 
tial, the averaging over waiting times (Eq- [8| affects D^. 
^ This was considered when fitting the experimental value 
Z)s = D{ts) (where ^ O.lt^) to the corresponding simula- 
tion value. Moreover, we took into account that the choice of 
Z)s affects and in turn 4 and = Z)(4). This procedure 
provides the friction coefficient close to the water-glass inter- 
face; (^0* ^ 1.4(^0 for P > (and = (^o for P = and the 
simulations). This implies an effective diffusion coefficient 
Dq /Do = (^o/<^o effective Brownian time jt^ = * /<^o 

in the experiments with P > 0, while in the simulations (and 
experiments with P = 0) Z)q = Z)o and = ^b- 

This procedure provides £(P) and corrects for radiation 
pressure effects, while fixing the diffusion coefficient D(t) at 
short and long times. Nevertheless, the intermediate, subdif- 
fusive behaviour is not affected and experimental and simula- 



tion results can be compared (Fig. 11 ). Agreement between 
experiments and simulations is found for D(t\ including in- 
termediate times. 

The diffusion coefficient at long times, Z)oo, has been 
linked to the free diffusion coefficient Z)o in the case of one- 
dimensional random potential energy landscapes with the val- 
ues of the potential following a Gaussian distribution with 





Fig. 12 Ratio of the long-time diffusion coefficient Doo and the 
short-time diffusion coefficient as a function of degree of 
roughness e in two different representations, as obtained from 
simulations (open squares) and experiments (filled circles, which are 
consistent with the simulation data by construction, Sec. |3.3| l. The 
lines represent (A) a spline interpolation and (B) a slope of two as 
expected for the case of a one-dimensional random potential energy 
landscape. 



width e:lSl 



(11) 



To our knowledge, theoretical predictions for higher dimen- 
sions are not available. Our simulation data (and, due to the 
procedure to link P io £, also the experimental data) sug- 
gest a similar dependence at small values of £ (Fig. [12^). 
For large e, however, we observe a significantly reduced de- 
crease of Doo /Do, i.e. a reduced slope in — In (Z)oo(e) /Do) and 
thus an exponent below 2 in equation ( pT) , maybe even a sat- 
uration. Therefore, for large £, the dynamics appear not to 
slow down as strongly as in a one-dimensional system.l^SIEIini 
This might be due to the possibility to avoid large barriers in 
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Fig. 13 Normalized diffusion coefficient D{t)/DQ sls sl function of 
delay time t for particles in a one-dimensional (lines) and 
two-dimensional (symbols) random potential energy landscape with 
different standard deviations e (as indicated) as observed in 
simulations. 



two dimenisons, which becomes increasingly attractive as £ 
grows. It might, however, also result from the fact that, with 
increasing e, the system requires an increasingly long time to 
reach equilibrium and hence the average over waiting times is 
strongly weighted toward non-equilibrium states (which also 
affects Z)s). Furthermore, with increasing £ it becomes in- 
creasingly difficult to determine Z)oo, i.e. the plateau of D{t) at 
large times. 



3.4 Comparison to particle dynamics in one-dimensional 
random and periodic potentials 

The particle dynamics in one-dimensional random potential 
energy landscapes also show three distinct regimes; diffu- 
sion at short and long times and subdiffusion at intermediate 
times (Fig. 13 )."^^^The dynamics are slower than in the two- 
dimensional case. In particular at long times, the dynamics 
are reduced and the long-time diffusive regime established at 
much longer times. In general, the characteristic times, for ex- 
ample ti_i and ta, are longer and show a stronger dependence 
on £. 

In a periodic sinusoidal potential, only one barrier height 
exists and thus the distribution of escape times is narrower. 
Thus, the dynamics at intermediate times are more suppressed 
with a smaller slope at the inflection point of the mean squared 
displacement, corresponding to a deeper minimum of the ex- 
ponent ld{t).^^ On the other hand, long-time diffusion is es- 
tablished earlier. 



4 Conclusions 

We investigated the dynamics of individual colloidal parti- 
cles in two-dimensional random potential energy landscapes, 
whose values follow a Gaussian distribution with a standard 
deviation e, which characterizes the degree of roughness of 
the potential. In the experiments, the potential was created us- 
ing an optical set-up and the roughness £ controlled via the 
laser power P. The relation £{P) could be established with 
the help of Monte Carlo simulations. The experimentally ob- 
served dynamics agree with our simulation results. Three dis- 
tinct regimes have been observed. At short times, the particles 
exhibit diffusive behaviour within their local minima, in which 
they remain until they cross a barrier, i.e. a saddle point, to a 
neighbouring minima. The wide distribution of barrier heights 
leads to a significant spread in residence times. This is re- 
flected in the mean squared displacement as a broad subdiffu- 
sive region with a relatively large slope at the inflection point 
at intermediate times. At long times, the hopping between 
minima resembles a random walk and diffusive dynamics are 
recovered. The long-time diffusion coefficient decreases with 
increasing degree of roughness £. For large £ this decrease 
is less pronounced than in one-dimensional potential energy 
landscapes. ESI This is attributed to the possibility to bypass 
large barriers in two-dimensions, which is particularly utilized 
if large barriers are prevalent, i.e. for large £. 

The system presented here can also serve as a well- 
controlled, tunable and easily observable model for other sys- 
tems, which either explore space or configuration space, i.e. a 
potential energy landscape. These systems include crowded 
systems, such as concentrated colloidal suspensions, super- 
cooled liquids, glasses, 35 living cells, but 
also complex potential e nergy landsc apes, such as those sug- 
gested in protein MdingJ^^^^^^^ 
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